arXiv:1509.05564v3 [hep-ph] 27 Jan 2016 


Asymmetric dark matter and effective number of neutrinos 

Teruyuki Kitabayashi* and Yoshihiro Kurosawa^ 

Department of Physics, Tokai University, 4^-1-! Kitakaname, Hiratsuka, Kanagawa, 259-1292, Japan 

(Dated: January 28, 2016) 

We study the effect of the MeV-scale asymmetric dark matter annihilation on the effective number 
of neutrinos at the epoch of the big bang nucleosynthesis. If the asymmetric dark matter x 
couples more strongly to the neutrinos v than to the photons 7 and electrons e~, <C or 

3> the lower mass limit on the asymmetric dark matter is about 18 MeV for Ned — 3.0. 
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I. INTRODUCTION 

One of the big motivations for considering the asym¬ 
metric dark matter (ADM) scenario [1] is that the abun¬ 
dance of baryon Df, and dark matter (DM) Dqm is ob¬ 
served to be close to each other Ddm ^ 5Df,. In an ADM 
model [2], there is a common mechanism that might give 
rise to the baryon asymmetry as well as the ADM asym¬ 
metry and we obtain DADM/Hb — iVx/ 9 b)iiTix/mp), 
where Dadm, Vh, ''^p ADM abun¬ 

dance, the baryon asymmetry, the ADM asymmetry, the 
ADM mass and the proton mass, respectively. From this 
relation, the natural scale for ADM is around 5 GeV, 
however, mass of the ADM can be as low as a few keV in 
some models [3, 4]. Moreover, the light non-asymmetric 
DM (we call it symmetric DM or simply DM), such as 
MeV to GeV DM, remains an elusive blind spot in the 
current underground searches [5]. 

If there are light ADM particles, the energy density 
of radiation at big bang nucleosynthesis (BBN) epoch 
could be changed. It is not necessary that the particle be 
a dark matter candidate to influence BBN. It only need 
be a thermal relic, e.g., a particle that was in thermal 
equilibrium with the standard model particles present 
when the temperature was compatible with the mass of 
the relic particle. 

These extra particles contribute to the unknown radi¬ 
ation content of the universe [6-8]. The energy density 
of the relativistic particles at BBN (as well as at time of 
cosmic microwave background (GMB) photons released) 
is usually expressed in terms of the effective number of 
neutrinos Ngg. For the standard cosmology, Nes = 3 
is expected. Including the effect of slight reheating of 
the neutrinos from early e“'"e“ annihilation, we obtain 
Veff = 3.046 [9]. 

The relation between the symmetric DM at BBN and 
the effective number of neutrinos is extensively studied 
in the literature [10-23]. The extra particles at BBN 
contribute to the effective number of neutrinos in the 
following two cases: 
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1. Direct contribution case: If an extra particle is light 
enough such as eV-scale sterile neutrino, it is re¬ 
garded as one of the radiation components. This 
light particle has contributed directly to the effec¬ 
tive number of neutrinos as the so-called dark ra¬ 
diation [14, 15]. 

2. Indirect contribution case: Although the extra par¬ 
ticles are not light enough to contribute directly to 
the effective number of neutrinos, its annihilation 
heats other particles via entropy transfer. Gonse- 
quently, these extra particles contribute to the ef¬ 
fective number of neutrinos indirectly even in the 
absence of dark radiation. Either increase or de¬ 
crease of the effective number of neutrinos occurs 
as follows: 

(i) If an extra particle couples more strongly to the 

neutrinos v than to the photons 7 and electrons 
e“, its late time annihilation heats 

the neutrinos more than the photons. Ultimately, 
this type of extra particle yields an excess of the 
effective number of neutrinos, N^g > 3 [10, 11, 18, 
19, 21]. 

(ii) On the contrary, if an extra particle couples 
more strongly to the electrons and photons than 
to the neutrinos, Tx-y,xe ^x^y if® annihilation 
heats the electron-photon plasma relative to the 
neutrino background, leading to a reduction in the 
effective number of neutrinos below the standard 
model value, Ngg < 3 [10, 11, 13, 18, 19, 21]. 

The relic density of ADM in the original models is set 
by the asymmetry around the time of baryogenesis to 
obtain the single explanation for both baryon and dark 
matter densities [24-27]. 

Some scenarios to generate primordial dark and baryon 
asymmetry in the ADM models are proposed such as 
lepton-number and/or baryon-number violating decay, 
the Affleck-Dine mechanisms, via phase transition at 
electroweak baryogengesis, scenarios with dark gauge 
group and messengers between dark and visible sector 
(see Ref.[l] and references therein). For some of these sce¬ 
narios, since both DM and anti-DM particles may popu¬ 
late the thermal bath in the early universe, the relic num¬ 
ber density of ADM is determined not only by the initial 
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number asymmetry but also the annihilation cross sec¬ 
tion [28-30]. This asymmetric WIMP framework can ac¬ 
commodate a wide range of dark matter masses and anni¬ 
hilation cross sections. For example, Graesser, et al. ap¬ 
plied the asymmetric WIMP scenarios to lepton-number 
violating and baryon-number violating ADM models [28]. 
Another example is given by Lin, et ah, they reported the 
model independent constraints on the quantities of light 
ADM with mass ~ 1 MeV - 10 GeV in the asymmetric 
WIMP framework [29, 30]. 

The generation mechanism of initial asymmetry, along 
with ongoing self-annihilation, would be important sub¬ 
ject in the study of ADM. We keep the similar mechanism 
which is proposed in Ref. [28] in mind and we put aside 
discussion of the particular mechanism to generate the 
initial asymmetry. We would like to show a model inde¬ 
pendent analysis in this asymmetric WIMP paradigm. 

In the ADM scenarios, the ADM particle x its 
antiparticle x have nonzero chemical potentials and the 
particle x is not self-conjugate of the antiparticle y. To 
study the dependence of the light ADM on the effective 
number neutrinos, we should take care of the chemical 
potential of ADM. 

We have a little knowledge about the dependence of 
the extra particle with nonzero chemical potential on the 
effective number of neutrinos. For the direct contribu¬ 
tion case, a constraint on the chemical potential for a 
fermionic light DM has been deduced from BBN calcu¬ 
lation by Boeckel and S.-Bielich [31]. The correlation 
between the effective number of neutrinos and the cos¬ 
mological parameters with light DM particles (including 
ADM) have been studied by Blennow et al. [32]. These 
two papers give us many interesting results related to 
the ADM, however, the dependence of the ADM number 
asymmetry on the effective number of neutrinos is not 
clear yet. For the indirect contribution case, there was 
no study of the relationship between the ADM and the 
effective number of neutrinos. 

In this paper, the known two methods to estimate the 
effective number of neutrinos with light extra symmetric 
DM by Boeckel and S.-Bielich in the direct contribution 
case [31] and by Steigman in the indirect contribution 
case [18] are slightly extended to the light ADM in a 
straightforward way. By using the extended methods, 
we discuss some constraints on a MeV-scale ADM with 
the effective number of neutrinos in both of direct and 
indirect contribution cases. Significant lower limit on the 
ADM mass is obtained in the indirect contribution case. 
The upper limit on the ADM number asymmetry is also 
obtained in both cases, however, we observe that this 
upper limit is not strongly constrained by the relic abun¬ 
dance consideration. 

This paper is organized as follows. In Sec. II, we review 
the basic picture of the ADM and the effective number of 
neutrinos. This review does not include any new findings. 
We would like to present a brief review of the ADM, a 
detailed analysis of the calculation of the relic abundance 
of a thermal dark matter candidate that presented in 


many previously published papers. We also show our 
notations in this section. 

The new results are reported in Secs. Ill and IV. In 
Sec. Ill, we show the method to obtain a constraint on 
the chemical potential of ADM in the direct contribution 
case which is developed by Boeckel and S.-Bielich [31]. 
We use their method to obtain the limit on the ADM 
number asymmetry with the effective number of neutri¬ 
nos. This is the first result in this paper. In Sec. IV, we 
show the useful method which is developed by Steigman 
[18] to estimate the effective number of neutrinos with 
symmetric DM at BBN in the indirect contribution case. 
Then, we extend this method to include ADM. This ex¬ 
tension is second and main result in this paper. In the 
same section, the constraints on the mass, number asym¬ 
metry and cross sections for ADM are obtained numeri¬ 
cally. The actual calculations are important complement 
to any previously reported results in the literature. These 
complementary results are the third result in this paper. 
Finally, Sec. V is devoted to a summary. 


II. ASYMMETRIC DARK MATTER AND 
EFFECTIVE NUMBER OF NEUTRINOS 


A. Relic abundance 


Relic abundance of the ADM has been studied [28- 
30, 33-38] based on the methods for the symmetric DM 
[33, 39-42]. We assume that, at the moment close to the 
ADM decoupling epoch, the only reactions that change 
the number of x and x are annihilations and pair creation 
of XX ^ f f (there is no self-annihilation and creation 
such as XX // and xx ^ ff [35]). With this assump¬ 
tion in mind, the relic density of y and y is determined 
by solving the following Boltzmann equations 
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where Uy and Uy denote the number density of y and 
y, respectively. Both y and y may populate the thermal 
bath (equilibrium) in the early universe. The equilibrium 
densities and in the presence of asymmetry dif¬ 
fer by the chemical potential Hy. The detail of the an¬ 
nihilation process is included by the thermally averaged 
annihilation cross section {ayyv). The Hubble expansion 
rate is calculated as 


H = 

Mpi V 90’ 


( 2 ) 


during the radiation dominated epoch where Mpi = 
2.4 X 10^® GeV, and T denote the reduced Planck 
mass, the effective relativistic degrees of freedom for the 
energy density and the temperature of the thermal bath 
(temperature of the photons), respectively. The effective 
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relativistic degrees of freedom is defined as 

4 


9*{T) = X! 




T 


7 


E 




2 =fermions 


, ( 3 ) 


where Ti and gt are the temperature and the number 
of internal degrees of freedom of species respectively 
[43, 44], 

We use the standard definitions 




X = 


T ’ 






where 


- ^ t3 

* “ 45 9*sT , 


is the entropy density and 


9*s(T) = 


9^ 


i=bosons 


E 9 ^ 

■i^fermions 


( 4 ) 


( 5 ) 


, (6) 


is the effective relativistic degrees of freedom for the en¬ 
tropy density for relativistic particles, which turns out to 
be 

g«(T)= ^ 9^F- + l Y. 9^Ft, (7) 

2 =bosons z=fermions 


in a more general case, where 
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( 8 ) 


is a function of the particle mass that change smoothly 
from = 1 when the particle is ultrarelativistic (xi = 
rrii/T <C 1) to = 0 when it becomes nonrelativistic 
(x, » 1) [32], 

We assume that the universe expands adiabatically and 
that g* as well as are treated as a constant during 
the XX annihilation period. In terms of and x, 

the Boltzmann equations become 


dY^ _ dY^ 
dx dx 


27r2 g*m^ / 

H 45 J’ 

(9) 


and we obtain c?(I^ — Y^)/dx = 0. Thus the ADM num¬ 
ber asymmetry, namely, the net comoving densities, 

e = r^-y^, ( 10 ) 

is constant. Because the number asymmetry e is con¬ 
served, the ADM asymmetry in equilibrium remains 

at any time. Using the definition of the number asym¬ 
metry e in Eq.(lO), we obtain the final form of the Boltz¬ 
mann equations 



{a^X^) 27r2 g,ml 

(YYeY^-YYY^^), 

= 

dx 

H 45 a;4 


(o’xx'v) 271"^ 9*m^ 

(Y^ + eEx-Y^Y^). 

dx 

H 45 a;4 


( 11 ) 


We express the present relic abundance of the particle 
X in the terms of the density parameter times the 
scale factor for the Hubble expansion rate h = 0.673: 

^ oo) ^^ 2 ) 

/^crit /^crit 

where sq = 2.89 x 10^ cm“3 and pcrit = ^^/{SttG) = 
1.05ft-^ X 10“® GeV cm“3 are the present entropy density 
and the present critical density, respectively. We use the 
relation of the Hubble expansion rate Hq = 100ft- km 
s“^ Mpc“^ with the Newtonian gravitational constant 
G = 6.67 X 10“^^ m3 kg“^ s“^ [45]. The present relic 
abundance of the ADM is calculated to be 

UADMft^ = D^ft^ -|- 

= 2.75 X 10«^yADM(x ^ oo), (13) 
GeV 

where 


UADM(a;) = Yxi^) + Yx{x). (14) 

The observed energy density of the cold dark matter 
component in ACDM model by the Planck Collaboration 
is DDMft^ = 0.1188 ± 0.0010 (68% C.L.) [46]. 


B. Asymmetry and chemical potential 

The distribution function of particle species i is given 

by 

= e(is,-E/U + 0 ^’ 

where gi, Ei, fj,i and Ti denote number of internal de¬ 
grees of freedom, energy, chemical potential and temper¬ 
ature of particle i, respectively. The discrete parame¬ 
ter Qi takes only the following three values: 0^ = -1-1 
and 0i = — 1 correspond the Fermi-Dirac distribution 
and the Bose-Einstein distribution while 0i = 0 corre¬ 
sponds the Maxwell-Boltzmann distribution. In this pa¬ 
per, we consider a fermionic Dirac type ADM y and we 
take 0;^- = -1-1. 

The number density n^, energy density pi, pressure Pi 
and entropy density Si of particle species i with mass rrii 
in the isotropic universe are obtained as follows [43, 47]: 

1 

n* = ^ y E{E^ - mjfFf.dE, 

1 poo 

F = ^ j - mjYFhdE, 

1 

PiE Pi - ytTit 

■ 

i 

If the degeneracy of the particle i is small {p,i <C 1)), 
the Fermi-Dirac distribution can be well approximated 


(16) 

(17) 

(18) 
(19) 
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by the Maxwell-Boltzmann distribution. In the approxi¬ 
mation of Maxwell-Boltzmann statistics, the equilibrium 
number density of the ADM particle x becomes 


n 


EQ 
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= 9x 


( mxTx 
V 27r 


3l2 





for the nonrelativistic limit, ^ ^x ( 2 ^ ^ 1)) ^^id 


for the nonrelativistic and 


45 ffx 

27r‘* ’ 


( 29 ) 


for the ultrarelativistic cases. 

’ We observed the ratio /i^/^x many equations. This 
( 20 <)atio is the so-called degeneracy parameter or the pseudo¬ 
chemical potential. In the remaining part of this paper, 
we use the following definition 


( 21 ) 

for the ultrarelativistic limit, nix ‘C (x <C 1) [47]. The 
number density of antiparticle x is also obtained with the 
fact that = —fJ^x in equilibrium. As a result, chemical 
potential drops out in the product 


(30) 

2 x 

to denote the chemical potential of the ADM particle x 
and call ^ chemical potential simply. 

C. Effective number of neutrinos 


y^EQy'EQ 
X X 


1 



( 22 ) 


for the nonrelativistic and 


The energy density of relativistic particles, i.e., radia¬ 
tion components in the early universe, prad is given by 

Prad = P 7 T T PDRj (31) 


yEQy-EQ 
X X 


TT^ \27r^J \g^s 


(23) 


for the ultrarelativistic cases. The absence of the chemi¬ 
cal potential in Y^^Y^^ simplifies the Boltzmann equa¬ 
tions in Eq.(ll) in the nonrelativistic and ultrarelativistic 
cases [29]. 

The chemical potential of the ADM in equilibrium is 
still surviving on the asymmetry 

^EQ ^ ^.EQ -yP = -^ - nf Q) . (24) 


where p-y = (7r^/30)p.yT.^ is the energy density of pho¬ 
tons and pt“ = Nf’^p^ = 7V^“(7/8)(7rV30)p^T4 is the 
energy density of standard-model massless neutrinos for 
neutrino families. By the simple estimation in the 
standard particle cosmology, we have = 3. Includ¬ 
ing the effect of slight reheating of the neutrinos from 
early e+e” annihilation, we obtain = 3.046, which 
is due to the small overlap of neutrino decoupling and 
e+e" annihilation [9]. 

The energy density of extra radiation, dark radiation, 
is generally parametrized through the number of extra 
effective neutrino species /S.Ny as follows [8] 


For the examples, we obtain 

gEQ ^ 

27r^ g^s \ 27rT;^ ) 


p^xl'^x _ p~^xl^-) 


Pdr = 


(32) 


, In terms of the following effective number neutrinos 
(25) iVeff = + AfV,, (33) 


for the nonrelativistic [36] and 

^ (26) 

for the ultrarelativistic cases. From these equations, we 
have [36] 


^ =ln 
Tv 


1 

2 



where 


45 Px 

27r2 



3/2 


(27) 


(28) 


the total energy density of radiation components is esti¬ 
mated to be 


/2rad — Pj T 


8 g-y 



Pj. 


(34) 


As we mentioned in the introduction, if there is an 
extra light particle as a dark radiation that particle in¬ 
creases the effective number of neutrinos (the direct con¬ 
tribution case) [7, 18]. On the other hand, although the 
extra particles are not light enough to be a dark radi¬ 
ation, its annihilation heats other particles via entropy 
conservation and changes the ratio of (T^/Tx) in Eq.(34). 
The extra particle contributes to the total energy den¬ 
sity of radiation components Prad- Consequently, the 
change of prad reflects to the enhance or the reduce of 



5 


the effective number of neutrinos A^eff via the relation of 
Prad = P 7 + NgffPi/ even in the absence of dark radiation 
(the indirect contribution case). 

The effective number of neutrino can be probed by its 
effect on the CMB and the outcome of BBN. The recent 
observational result on the effective number of neutrinos 
by the Planck Collaboration is Nee = 3.04 ± 0.18 ( 68 % 
C.L.) from CMB data [46]. On the other hand, we have 
Nee = 3.71;[];4^ from BBN data [ 6 ]. 

III. DIRECT CONTRIBUTION 
A. Methods 



The method to obtain a constraint on the chemical 
potential of fermionic DM with the effective number of 
neutrinos in the case of direct contribution is developed 
by Boeckel and S.-Bielich [31]. First, we review their 
method and then we slightly extend their approach to 
estimate the upper limit on the ADM number asymmetry 
with the effective number of neutrinos. 

Because the energy density of ADM cannot excess 
the energy density of dark radiation pjyn [31], we find the 
following constraint on the energy density of ADM 

^ J A[/^(p,t) + /^(p,t)]d3p< (35) 


where fi {p, t) is the distribution function of ADM {i = 
X, x). In general, the energy density of fermion increases 
with its chemical potential. We obtain the upper bound 
on the chemical potential of ADM, as a function of 
the excess of the effective number of neutrinos at 

BBN, e.g., = /(AiV,). 

We slightly extend this method to estimate the upper 
limit on the ADM number asymmetry with the effec¬ 
tive number of neutrinos. Because the number density of 
fermion increases with its chemical potential, the upper 
bound on the net ADM number density, = n^ — n^ is 

a function of the excess of effective number of neutrinos 
ANi, at BBN, e.g., = /(AAj,). We can calculate 

the upper bound of the ADM number asymmetry at BBN 


BBN ^ n^^-jAN,) 
^max 5 


(36) 


where s is the total entropy density s = ^^Si. Recall 
that the number asymmetry is conserved, the upper limit 
of the ADM number asymmetry at any time Cmax j at the 
epoch of equilibrium and at BBN e®; is the same 


^max 

Qg g — gEQ _ BBN 
cmax Cjjjax '^max 


B. Constraints 


FIG. 1: The dependence of the extra effective number of neu¬ 
trinos ANi, on the maximum ADM number asymmetry Cmax- 


and we have 


9x 


— + k^ + —e 

120 4^ 87 r 2 ^ 


< AN 


III 



(38) 


from the relation of < poR- The upper bound of the 
chemical potential p,^ (more precisely degeneracy param¬ 
eter ^) of the ultrarelativistic ADM is obtained [31]: 


eBBN 

Smax 





The number density of ultrarelativistic ADM is estimated 
as = ( 5 xT'^/ 67 r^)( 7 r^ -|- ^^) and the upper bound of 
the ADM number asymmetry is calculated to be 

_ BBN _ 9x^x r^BBN I /aBBN\31 
tmax ^max %'k‘^S ' V'^max / J ■ 

The ratio of the temperature of neutrinos to the tem¬ 
perature of ADM, Ty/T^, can be estimated by consider¬ 
ing conservation of the comoving entropy. At BBN, we 
have [31] 


Ws{Tbbn)) ’ ^ ^ 

where T^d is the ADM decoupling temperature and Tbbn 
is the BBN temperature around 1 MeV. The effective 
degrees of freedom at BBN for standard model particles 
is obtained as 5 *s(Tbbn) = ‘2+ (7/8)(2 • 2 -|- 2 • 3) = 10.75 
for photons, and three neutrinos. Including ADM, we 
have 


For the ultrarelativistic ADM with temperature T^, 
the energy density is obtained as 


Pxx - Px'^x ( + 7^^ + 772^"^ 


1 


(37) 


<?*s(Tbbn) = 10.75 + ^-9x{FI[ + F+), (42) 

where is given by Eq.( 8 ). We take F+ = = 1 and 

Px = 2 for the ultra-relativistic Dirac ADM, (jf»s = 14.25. 
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The effective degrees of freedom at ADM decoupling 
9*s{Txd) depends on the details of the model of the 
ADM [31]. If the ADM is decoupled when all standard 
model particles were present and in equilibrium, we ob¬ 
tain (7*s(7x(i) = 106.75. In this section, we treat g^s^T^d) 
as a free parameter. 

Figure 1 shows the dependence of the extra effec¬ 
tive number of neutrinos /S.Ni, on the maximum ADM 
number asymmetry Cmax with rrix = 1 MeV. In the 
case of Cmax ^ 0 . 01 , the dependence of on the 

Cmax is almost negligible, e.g., we obtain the plateau for 
Cmax "S 0.01 in each curves. The constraint on e with 
is obtained: 


en,ax 0.01, (43) 

for AA^ < 2.67 and g^siT^d) > 10.75. 

In the appropriate GeV-scale ADM models, small num¬ 
ber asymmetry e~ 10 “^^ is favorable to explain the ob¬ 
served dark matter density Ddm (for example, see [29]). 
For the MeV ADM, the small number asymmetry is still 
expected (we will consider the relic abundance in Sec. 
IVB). The result of Cmax — 0.01 is not a strong con¬ 
straint. The cosmological consequences of on some 
problems are important [31], however, we leave the dis¬ 
cussion of the number asymmetry of the ultrarelativistic 
ADM in the direct contribution case. 

Because the nonrelativistic particles role as matter 
components, the direct contribution of such particles to 
the energy density of radiation components should be 
negligible. Consequently, we can hardly probe the chem¬ 
ical potential as well as the asymmetry for nonrelativistic 
ADM via BBN [31]. 


IV. INDIRECT CONTRIBUTION 

A. Methods 

The method to obtain the constraint on the effective 
number of neutrinos with the symmetric MeV DM in the 
indirect contribution is developed in the literature [ 10 , 11 , 
13-15, 18, 19, 21, 31]. First, we review the useful method 
for the symmetric MeV DM developed by Steigman [18] 
and then we extend his method to study the ADM in a 
straightforward way. 

Symmetric DM case (i) ^ the 

symmetric MeV DM y couples more strongly to the neu¬ 
trinos ly than to the photons 7 and electrons e~, e.g. 
Tx'y,xe F^^, its late time annihilation heats the neu¬ 
trinos more than the photons while the annihilation of 
the pairs heats the photons (but not the decoupled 
neutrinos). The entropy in the comoving volume of pho¬ 
tons and pairs S^e = B?{s-y -l-Se) and of the neutrinos 
and the symmetric MeV DM S^x = + s^) are con¬ 

served individually where R denotes the scale factor [43]. 
The ratio of neutrino to photon temperatures after both 
of the and the symmetric MeV DM are obtained by 


consideration of entropy conservation as follows: 

g~f ge4'e 1 -|- 

(44) 

where gi denotes effective internal degrees of freedom of 
the particle i; gi = (7/8) • 2 • 2 = 7/2 for the Dirac 
fermions, gi = (7/8) •2-1 = 7/4 for the Majorana 
fermions. Although we consider the only fermionic Dirac 
type DM, we can also use gi = I for the scalar bosons 
and 3 for the vector bosons with the bosonic distribu¬ 
tion function (Eq.(15) with 0^ = —1). The function of (j) 
denotes the normalized entropy density 


4>a{x) 


s“t(0) ’ 


(45) 


where Sa®*(a;) is the net entropy density of particle species 
a with vanishing chemical potential [18]: 

s“‘(2;) = Sa{g.a = 0) -k Sa{g,a = 0) = X! ^ ■ 

J-i 

z=Q;,a 

(46) 

For pairs, (j)e is evaluated at a; = Xed = "mejTvdi while 
for the symmetric DM, (px is evaluated at x = Xxd = 
fnxjTvd where denotes the decoupling temperature 
of neutrinos. Because the difference between Majorana or 
Dirac nature of symmetric MeV DM is taken into account 
by the effective internal degrees of freedom we take 
= 1 in Eq.(15) through our calculations. 

With the appropriate assumptions of Ty^ = 2 MeV and 
(pe = 0.993 [18], we obtain 


Aeff = 3 


11 (Ty 

t; 


4/3 


= 3.018 1 -f 




21 


4/3 


(47) 

If the symmetric MeV DM particle is sufficiently mas¬ 
sive, this particle is regarded as a matter component and 
Pxd ~ 0. As a result, the usual result NeS — 3 is recov¬ 
ered. On the contrary, for the very light symmetric DM, 
we obtain Px = sTi^xd )/~ s"®*( 0 )/s"®*( 0 ) = 1 
and the excess of the effective number of neutrinos N^g > 
3 without the dark radiation. The effective number of 
neutrinos is a function of the mass of symmetric MeV 


DM. 


Symmetric DM case (ii) Vxi,xe ^ ^xiy' On the 
other hand, if the symmetric MeV DM couples more 
strongly to the electrons and photons than to the neu¬ 
trinos, Tx'y^xe ^ ^xiy’ annihilation heats the electron- 
photon plasma relative to the neutrino background. By 
the consideration of entropy conservation, we have 



_57_ 

g-y + gepe + 9 x4’x 


2 

2 -I- l^pe + gx4>x 


(48) 


and 


A^eff = 3 


11 fn 


4/3 


= 3 


11 


10.95 ‘IgxPx 


4/3 


( 49 ) 
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If the symmetric MeV DM particle is sufficiently massive, 
the usual result A^eff — 3 is recovered. On the contrary, 
an annihilation of more light symmetric MeV DM reduces 
the effective number of neutrinos to be N^ff < 3. 

ADM case: We extend the method of calculating NeS 
with symmetric MeV DM in a simple way. To take care 
of the chemical potential of the ADM for the entropy 
calculations, we use the following net entropy density of 
the ADM 


formula in Eq.(54) to calculate the number density iii, 
we estimate the following integral [49, 50] 

ni = — / Ti{TiX + mi)[{TiX + mi)‘^ 

(55) 

instead of Eq.(16), where x = {Ei — mi)/Ti and the dis¬ 
tribution function in Eq.(15) is replaced by 


f-O = »> + »>= V (50) 

instead of Eq.(46). The normalized entropy density of 
the ADM to be a function of not only x^d but also is 
given as follows: 


4'xi^xd^ t^x) 


sT(^Xd,t^x) 


(51) 


Other quantities in Eq.(44) and Eq.(48), such as gx,4’e, 
remain the same. 

The effective number of neutrinos with the MeV ADM 
is estimated in the same form of Eq.(47) and Eq.(49). If 
the MeV ADM y couples more strongly to the neutrinos 
v than to the photons 7 and electrons e~ {Xxi,xe ^x’')^ 
we obtain 

Aeff = 3.018 (^i + ^sM^xd,f^x) y\ (52) 

while if the MeV ADM couples more strongly to the elec¬ 
trons and photons than to the neutrinos (T^-^^^e ^ 
we have 

/ 11 

^\l0.95 + 2gx(j)xix^d,g.x)) ' 


These derivations of the effect of light ADM on the 
effective number of ultrarelativistic species at the epoch 
of BBN, Eqs.(52) and (53), are the main new findings in 
this paper. 

In the following numerical calculations, we vary the 
mass of the MeV ADM in the range of = 0.1 — 100 
MeV. In this case, both of the ultrarelativistic limit and 
the nonrelativistic limit are inappropriate. We estimate 
the thermodynamic quantities such as the number den¬ 
sity from the first principle (most fundamental formulae) 
in Eq.(16), Eq.(17), Eq.(18) and Eq.(19) numerically by 
the Gauss-Laguerre integration method: 


N 

x°‘e~^F{x)dx = ''^WjF{xj), (54) 

where F(x) and Wj denote a function of x and “Gauss- 
Laguerre weights,” respectively [48]. In our calculation, 
a = 0 and A = 50 are chosen. In order to apply the 
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g(TiX+mi-^ii)/Ti _j_ 0^ 


(56) 


Similarly, the energy density pi and the pressure Pi are 
calculated as 

Pi = J Ti{TiX + mif[{TiX + mif-m^]^^'^fidx, 


1 


P, = —l T,[{T,x + m,f-mlf/^hdx, (57) 


instead of Eq.(17) and Eq.(18), respectively. 

Before we show the results from the numerical calcula¬ 
tions, we would like to comment on the role of asymmetry 
at BBN. The effect of the annihilating thermal relic par¬ 
ticles to heat either the photons or the light particles. 
The bulk of this heat is created when the temperature is 
within a factor of a few of the relic particle mass T ~ rriy,. 
At these temperatures the number asymmetry between 
particles and antiparticles may be small without large 
chemical potential. At these temperatures the number 
of particle and antiparticle pairs could be much larger 
than the relic number when annihilations have ceased, at 
a temperature much lower than the particle mass. Thus, 
the effect of the number asymmetry of unobserved MeV 
scale particles is not strong at BBN (except the very light 
and degenerate (asymmetric) neutrinos. Eor example, see 
Refs.[51, 52]). 

Although the case of asymmetric dark matter is not 
significantly different from the case of symmetric dark 
matter without large chemical potential of ADM, we 
would like to announce that the explicit considerations 
of the relation between the chemical potential of ADM 
and effective number of neutrinos at BBN are first shown 
in this paper. Some numerical results are expected with¬ 
out actual calculations, however, the actual calculations 
are important complement to any previously reported re¬ 
sults in the literature. 


B. Constraints 

1. Lower mass limits: Eigure 2 shows the depen¬ 
dence of the effective number of neutrinos Aeff on the 
mass of the ADM with the various chemical poten¬ 
tials ^ at BBN (T = 2 MeV). The figure (a) shows the 
dependence in the case of ^xi,xe ^xi' while the figure 
(b) shows the dependence in the case of ^xi,xe ^ ^xi'- 
The value of Ngf f with the vanishing chemical poten¬ 
tial (^ = 0 ) is compatible with the previously reported 




m^[MeV] 



FIG. 2: The dependence of the effective number of neutrinos 
A'eff on the ADM mass m-)^ MeV with the various chemical 
potential ^ at BBN. The figure (a) shows the dependence 
in the case of <C while the figure (b) shows the 

dependence in the case of 


results for the symmetric MeV DM case [12, 18, 21, 53]. 
For example, Boehm, et al. obtained the lower bound 
on the mass of Dirac fermion DM from the CMB data 
(and BBN considerations) at 95% C.L. as nix > 7.3 
MeV for Neff = 3.301q;5j = 2.79 — 3.84 in the case 
of obtain the same lower limit of 

rux for Neff = 3.86 from Fig.2. We note that the value 
of Neff at CMB and at BBN may be different, e.g., a 
MeV scale particle with the vanishing chemical potential 
(mass < 10 MeV) may lead to the value of Neff at CMB 
formation which is lager than at BBN [12]. In this paper, 
we take Neff to be defined at BBN. 

The similar figure of Fig.2 has already reported by Nol- 
lett and Steigman for symmetric DM [19, 21] . The curves 
with nonvanishing chemical potential in Fig.2 are newly 
obtained in our study. The dependence of the chemi¬ 
cal potentials of ADM on the effective number of neu¬ 
trinos is shown explicitly for the first time in this pa¬ 
per. Compere with the symmetric DM case, the effec¬ 
tive number of neutrinos Veg increases with the increas¬ 
ing asymmetry (chemical potential ^) and with the de¬ 
creasing mass rUx if ADM particle mainly interacts to 


neutrinos (F^-^^^e 'C ^x^)- the contrary, Vgg de¬ 
creases with the increasing ^ and with the decreasing 
if ADM particle mainly interacts to photons and elec¬ 
trons Txi^)- Also, in the Nollett and Steigman 

papers for symmetric DM [19, 21], Vgg depends on the 
nature of the quantum statistics of thermal relic (i.e., 
fermion or boson). We can take the effective internal de¬ 
grees of freedom as gx = 1,7/4, 2,7/2 for a real scalar, 
Majorana fermion, complex scalar and Dirac fermion in 
Eqs(52) and (53). Thus, N^a depends on the difference of 
the nature of the thermal relic not only in the symmetric 
DM case but also in the ADM case. 

The lower mass limit is obtained with the upper bound 
or lower bound on the effective number of neutrinos. For 
example, we obtain rUx ^ = 18.1 if N’^^ = 3.0 

in the figure (a), while = 18.3 if = 

3.0 in the figure (b). Thus, the sets of for 

the fixed or is obtained. The lower bound 

on the ADM mass is the case when the asymmetry and 
chemical potential are small. The bound gets stronger 
as asymmetry e and chemical potential ^ grow. It shows 
the smooth transition from asymmetric to the symmetric 
WIMP limit. 

The sets of {(^,m™™)} is shown in Fig.3. This figure 
shows the dependence of the lower mass limit of the ADM 
on the chemical potential ^ with the various bound 
of the effective number of neutrinos Vgg at BBN. The 
figure (a) shows the while the figure 

(b) shows the ^xi,x<i ^ case. The ADM mass m™'" 
increases with the decreasing upper bound of the effective 
number of neutrinos N^^ in the Tx'y,x<^ case. On 

the contrary, m™'" increases with the increasing the lower 
bound of in the Tx-y,xe ^ case. Moreover m™'" 
increases with the increasing chemical potential ^ in both 
cases. 

From Fig.3, the following constraint on with N^s 
is obtained: 

8.1-18.1 MeV (V/l|“ = 3.0 - 3.7) 

8.2 - 18.3 MeV {N^'^ = 2.5 - 3.0) ^ ’ 

for Tx-y^xe ^xi' (cipper) and ^ (lower) 

cases, more concretely, 

18.1 MeV (Tx^,xe « fVgT" = 3-0) 

18.3 MeV {rxj,xe » N^^ = 3.0) ^ ’ 

with the standard value of the effective number of neu¬ 
trinos. 

2. Maximum number asymmetry: Once the sets 
of {(^jW™'")} is obtained, we can estimate the ADM 
number asymmetry at (^,m™™): 

e(C, m““) = i [n^(e, - nx{^, ")] . (60) 

The total entropy density s is calculated as s = Ssm+'^x*** 
where ~ s"®*-|-Sg®*-|-3s])®* denotes the sum of the net 
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FIG. 3: The dependence of the lower mass limit of the ADM 
on the chemical potential ^ with the various effective 
number of neutrinos Nes at BBN (T = 2 MeV). The figure 
(a) shows the <C T^^^ case while the figure (b) shows 

the f S> case. 


entropy densities of the standard model particles at BBN 
and denotes the net entropy density of the ADM 
y. In the numerical calculation of the entropy density 
[see Eq.(19) with Eq.(55) and Eq.(57)], the masses of 
all particles except neutrinos are taken from the parti¬ 
cle data group [45]. Although, the recent results of the 
Planck experiment give us a constraint on the sum of the 
light neutrino masses as < 0.17 eV (95% C.L.) 

[46], we assume that all neutrino masses are the same 
as TOy = = rriv^ = = 1 eV in our calculation 

for the sake of simplicity. The unknown tiny masses of 
the neutrinos are irrelevant in our study. Moreover, we 
neglect the chemical potential of all standard model par¬ 
ticles. 

Because the ADM number asymmetry e decreases with 
the increasing mass (we will see below), we obtain the 
upper limit of the ADM number asymmetry 

e„.ax = e(e,<“), ( 61 ) 

and the set of {(^, Cmax)} for the fixed or 

Figure 4 shows the ADM number asymmetry. Figure 
(a) describes the dependence of the ADM number asym- 



[MeV] 




FIG. 4: ADM number asymmetry. Figure (a): the depen¬ 
dence of the ADM number asymmetry e on the mass 
with fixed chemical potential Figure (b) and (c): the de¬ 
pendence of the upper limit of the ADM number asymmetry 
Emax on the chemical potential ^ with the various bound of 
the effective number of neutrinos or at BBN. The 

figure (b) shows the case while the figure (c) 

shows the ^x-y,xe ^x<' case. 


metry e on the mass with fixed chemical potential 
Figures (b) and (c) describe the dependence of the 
upper limit of the ADM number asymmetry Cmax on the 
chemical potential ^ with the various bound of the ef¬ 
fective number of neutrinos or at BBN. The 

figure (b) shows the case while the figure 
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(c) shows the ^ case. In the figure (a), we 

see that the ADM number asymmetry e decreases with 
the increasing mass m-^. The figure (b) shows that the 
upper bound of the ADM number asymmetry Cmax in¬ 
creases with the increasing effective number of neutrinos 
in the ^ F^^^ case. On the contrary, the 

figure (c) shows that Cmax increases with the decreasing 
in the ^ T^i, case. 

From Fig.4, the following constraint on Cmax with iVegf 
is obtained: 


^max — 

f 1.35 X 10-^ - 7.20 X 10-® = 3.0 - 3.7) 

11.23 X 10-^ - 6.95 X 10-® = 2.5 - 3.0) 

for ^ > 0.0025, (62) 

for <C Ty^ (uppor) and ^y-y,ye ^ (lower) 

cases, more concretely, 


€max 


1.35 X 10-7 (F^^.^e « ry,,N^s^- = 3.0) 

1.23x10-7 (Fx^.xe»rx.,lvmm^3 0 )l ’ 


with the standard value of the effective number of neu¬ 
trinos. 

We note that, as shown in figure (b) and figure (c), Cmax 
increases with the increasing chemical potential It is 
naturally expected because the origin of the ADM num¬ 
ber asymmetry is the nonvanishing chemical potentials 
of X Slid X- Although there is no indication of the upper 
limit of the ADM number asymmetry Cmax for ^ = 0 in 
Fig.4, Cmax = 0 is obtained correctly for ^ = 0. 

As we show below, in consideration of the relic abun¬ 
dance of the MeV ADM, the obtained upper limit on the 
ADM number asymmetry in Eq.(62) and Eq.(63) is not 
a strong constraint, if the origin of the observed 511 keV 
gamma-ray is the annihilating light MeV ADM. 

3. Relic abundance: The relic abundance of the 
ADM depends on not only the ADM number 

asymmetry e but also the thermally averaged annihilation 
cross section {cTyyv) [see Eq.(ll)]. The averaged cross 
section is obtained by [42] 

= smy^A-K.) <®'‘> 

pOO 

X / Crxx(s)('S - 4TO^)\/sA:i(Vs/T:()ds, 

J 4m| 

where Uyy, Ki, s and x denote the annihilation cross 
section, the modified Bessel function of order i, one of 
the Mandelstam variable (not an entropy density) and 
X = mylTy, respectively. Practically, the averaged cross 
section can be expanded in power of the relative velocity 
of incoming particles v. The standard approximation of 
the averaged cross section is [33, 35]: 

{ayyv) ~ a + bx~^+0{x~'^). (65) 

If the s-wave annihilation is dominant, we can take a ^ 0 
and 6 = 0. On the other hand, if the p-wave annihilation 





[MeV] 


FIG. 5: The dependence of the relic abundance of the MeV 
ADM Dadm^^ on the mass niy in the s-wave dominant case. 
The figure (a) shows the dependence with the various num¬ 
ber asymmetry e, while the figure (b) shows the dependence 
with the various cross section {cfyyv) ~ a -|- bx~^. The figure 
(c) shows the value of parameter “a” required to satisfy the 
observed relic abundance of DM, Dadm 6^ = 0.12 and lower 
mass limit for the N^ff = 3.0. The vertical line in figure (c) 
shows the lower mass limit. 


is dominant, we can put o = 0 and 6 7 ^ 0. In this study, 
the coefficients a and b are taken as free parameters to 
perform a model independent analysis. 

Eigures 5 and 6 show the dependence of the relic abun¬ 
dance of the MeV ADM on the mass niy in the 

s-wave dominant case (Fig.5) and in the p-wave domi- 
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mj( [MeV] 

FIG. 6: Same as Fig.5 but p-wave dominant case. 


nant case (Fig.6), respectively. In Fig.5, the figure (a) 
shows the dependence with the various number asymme¬ 
try e, while the figure (b) shows the dependence with 
the various cross section {o'xx'^) — a + bx~^. The figure 
(c) shows the value of parameter “a” required to satisfy 
the observed relic abundance of DM, Dadm^^ = 0.12 
and lower mass limit for the N^ff = 3.0. The vertical 
line in figure (c) shows the lower mass limit. The relic 
abundance of the MeV ADM Dadm^^ decreases with the 
decreasing ADM number asymmetry e as shown in the 
figure (a). The mass dependence on the Dadm^^ is sig¬ 
nificant for the very light uix ^ 5 MeV and the highly 
asymmetric e>lxl0“^ case. In other words, we can 


ignore the mass dependence on the relic abundance for 
the case of e < I x 10“®. Indeed, for e < 1 x 10“®, the 
s-wave dominant cross section is almost degenerate for 
= 0.12 as shown in figure (c). The behavior of 
the relic abundance in the p-wave dominant case (Fig.6) 
is similar to the s-wave dominant case (Fig.5). 

The cosmological constraints on the annihilation cross 
section of the symmetric DM are extensively studied in 
the literature, such as the constraints from the galactic 
511 keV gamma-ray line [54-64], from the galactic 3.5 
keV x-ray line [65-67], from other cosmic rays [68-73], 
from the CMB [74-77] and from the BBN consideration 
[78]. 

For the symmetric MeV DM, the significant constraints 
on the annihilation cross section are obtained from the 
galactic 511 keV gamma-ray line observation. The light 
particles y {rrix ^ 1 — 100 MeV) annihilating into e+e" 
pairs in the galactic bulge may be the source of the ob¬ 
served gamma-ray [54]. According to this scenario, if the 
s-wave annihilation is dominant, rux 100 MeV is re¬ 
quired to obtain the observed cosmic ray flux without the 
final state radiation/bremsstrahlung, the so-called inter¬ 
nal bremsstrahlung (IB), effects. Including IB effects, we 
may expect that the upper bound of the symmetric dark 
matter particles is around 3 MeV [79] or 7.5 MeV [80]. 
On the contrary, if the cross section is p-wave dominant, 
the lighter DM (rrix ~ 1 MeV) is possible for b ~ 10 — 200 
pb {b - 2.6 X 10-® - 5.1 X lO”'^ GeV^) [54, 58]. 

If the s-wave annihilation is dominant, the upper 
bound of symmetric dark matter rux 3 or 7.5 MeV 
is at almost margin of our result in Eq.(58) for N^ff = 
2.5 — 3.7. Moreover, for N^ff = 3.0, ^ 3 or 7.5 

MeV is not consistent with our result to™™ ~ 18 MeV 
in Eq.(59). Although, we may expect that the constraint 
on the symmetric MeV DM is also appropriate for the 
MeV ADM approximately, there is a little study for the 
cosmological constraints on the annihilation cross section 
of the ADM [81]. The upper mass limit of ADM may be 
a few MeV, however, we take a conservative upper limit 
of the mass of ADM as a few 10 MeV in this paper. 

Eigure 5 shows that the s-wave dominant annihilation 
with a ~ 8 X 10“® GeV-2, 6 = 0 and e ~ 1.0 x 10“® is 
consistent with the observed energy density of the cold 
dark matter component ~ 0.1. For the p-wave 

dominant annihilation, a = 0, 6 ~ 3 x 10“^ GeV“^ and 
e ~ 1.0 X 10“® is consistent with data as shown in Fig.6. 
If the origin of the observed 511 keV gamma-ray is the 
annihilating MeV ADM and/or all of the cold dark mat¬ 
ter is made of the MeV ADM Dadm^^ = Ddm^^, the 
ADM number asymmetry e should be less than the ob¬ 
tained upper limit Cmax in Eq.(62) and Eq.(63) at least 
under 10“^. If the effective number of neutrinos is just 
the standard value = 3.0, the ADM number asym¬ 
metry should be e ~ Cmax x 10“®. As a consequence, the 
upper limit on the ADM number asymmetry in Eq.(62) 
and Eq.(63) is not strong constraint. 

We comment that the constraints in Eq.(62) and 
Eq.(63) may be more strict in the combination of the 
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following cases: (a) if the chemical potential of the ADM 
is tiny ^ 0.0025, (b) if the origin of the observed 511 

keV gamma-ray is not the annihilating MeV ADM and 
(c) the MeV ADM is a part of the cold dark matter 

V. SUMMARY 

We have extended the known two methods by Boeckel 
and S.-Bielich [31] as well as by Steigman [18] to obtain 
the constraints on the MeV asymmetric dark matter with 
the effective number of neutrinos Ngg at big bang nucle¬ 
osynthesis. 

If an extra particle is light enough, this light particle 
has contributed directly to the effective number of neu¬ 
trinos as the so-called dark radiation (direct contribution 
case). From the requirement of < PdR) we have ob¬ 
tained the upper limit on the asymmetric dark matter 
number asymmetry Cmax ~ 0.01 for < 2.67. 

Although the extra particles are not light enough to 
contribute directly to the effective number of neutri¬ 
nos, its annihilation yields either increase or decrease 
the effective number of neutrinos (indirect contribution 
case). If the MeV asymmetric dark matter couples more 
strongly to the neutrinos than to the photons and elec¬ 
trons (r^.y_^e if MeV asymmetric dark 

matter couples more strongly to the electrons and pho¬ 
tons than to the neutrinos (Xxi,xe constraint 

on TO™'" with Veff is obtained in Eq.(58). For example, 

TO™™ ~ 18 MeV, (66) 

if the effective number of neutrinos is just the standard 


value Veff = 3.0 as shown in Eq.(59). The constraint 
on e with Ves is also obtained in Eq.(62). Eor example 
Emax — lO”"^ for Aeff = 3.0 as shown in Eq.(63). 

If the origin of the observed 511 keV gamma-ray is the 
annihilating MeV asymmetric dark matter and/or the 
all of cold dark matter is made of the MeV asymmetric 
dark matter, e ~ 1.0 x 10“® is consistent with data. The 
constraint on e is not a strong constraint in both of direct 
(fmax 0.01, for AVi, < 2.67) and indirect contribution 
cases (cmax — 10“^ for Aeff = 3.0). 

On the other hand, the lower limit of the asymmetric 
dark matter mass to™'" ~ 18 MeV for ~ 3.0 in the 
indirect contribution case is strict constraint. From the 
galactic 511 keV gamma-ray line observation, the sym¬ 
metric dark matter mass may be less than about a few 
10 MeV (or a few MeV). We can expect that the range of 
the asymmetric (as well as symmetric) MeV dark matter 
mass is so narrow to satisfy ~ 18 MeV in the case of 
f^x 7 .xe ^ ^x^' f^x7,xe Constraint may be 

useful to check or construct the MeV asymmetric dark 
matter models. 

The role of the number asymmetry of unobserved MeV 
scale particles is not strong at BBN. The case of ADM is 
not significantly different from the case of symmetric DM 
without large chemical potential; however, the explicit 
considerations of the relation between the chemical po¬ 
tential of ADM and effective number of neutrinos at BBN 
are shown for the first time. Our actual calculations are 
important complement to any previously reported results 
in the literature. The method to discuss the dependence 
of the chemical potential of ADM on the effective num¬ 
ber of neutrinos give a little useful step toward solving 
puzzles in the dark matter problems. 
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